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A isobar components in the nuclear many-body wave function are investi- 
gated for the deuteron, light nuclei ( 16 0) and infinite nuclear matter within 
the framework of the coupled-cluster theory. The predictions derived for var- 
ious realistic models of the baryon-baryon interaction are compared to each 
other. These include local (V28) and non-local meson exchange potentials 
' (Bonn2ooo) but also a model recently derived by the Salamanca group ac- 

counting for quark degrees of freedom. The characteristic differences which 
are obtained for the iVA and AA correlation functions are related to the 
approximation made in deriving the matrix elements for the baryon-baryon 
[ interaction. 

o 

5 1 I. INTRODUCTION 

o _ 

It has always been one of the basic challenges of theoretical nuclear physics to derive the basic 
properties of nuclei from realistic models of the nucleon-nucleon (NN) interactions. Very 
sophisticated approximation schemes have been developed to solve the many-body problem 
of strongly interacting nucleons. As examples we mention Brueckner hole-line expansion 
I]], the coupled-cluster or "exponential S" approach 0, the self-consistent evaluation of 
Greens functions [[J, variational approaches using correlated basis functions [Q and recent 
developments employing quantum Monte-Carlo techniques [[J. (see e.g. the recent review 
d | @ on these many-body approaches). 

A major ingredient for such investigations is the definition of a realistic NN interaction. 
Here we define a realistic NN interaction to be a model for the nucleon nucleon interaction 
in which the parameters have been adjusted to obtain a very good fit of the experimental NN 
scattering data at energies up to the threshold for n production. As examples for modern 
realistic interactions we refer to refs. f7[ [| . 

The central assumption of such investigations is that the nucleons can be considered as 
inert particles interacting by two-body forces, which are identical in the vacuum and in 
the nuclear medium. It is of course well known that nucleons are not elementary particles. 
At short distances, nucleons will polarize each other, which will lead to virtual excitations. 
These sub-nucleonic degrees of freedom are effectively taken into account in the realistic NN 
interaction as the parameters are adjusted to the data. 

As an example we mention the processes of the type that two nucleons interact with each 
other exchanging a n or p meson leading to a AA or A A state. After the exchange of a second 
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meson the two baryons may be returned into a NN state such that this process contributes 
to the amplitude for elastic NN scattering. A theoretical approach which accounts for 
interacting nucleons only parameterizes this process with intermediate isobar excitations by 
adjusting e.g. the coupling constant of the a meson in a One-Boson-Exchange model for the 
NN interaction [TI| . If, however, the isobar configurations are included explicitely one finds 



that a large part of the medium range attraction in the NN interaction originates from the 
coupling to the NA and AA configurations. 

Both kinds of models, with and without the explicit consideration of isobar configurations, 
yield essentially the interaction for two nucleons in the vacuum, since they are adjusted to 
reproduce the NN phase shifts. Remarkable differences show up, however, if the differ- 
ent kinds of interactions are considered for the interaction of two nucleons inside a nuclear 
medium. While the exchange of the a meson is essentially the same in the vacuum and in 
the nuclear medium, the formation of virtual iVA and AA configurations is suppressed in 
nuclear matter (due to Pauli effects and binding of the nucleons) as compared to the vacuum. 
This leads to a loss of attraction for the effective interaction of two nucleons in the medium, 
which results in a smaller calculated binding energy. This effect has been investigated by 
various groups using the Brueckner-Hartree-Fock approximation, which means that corre- 
lations between nucleons are taken into account by means of the effective interaction, the 



G-matrix [I I -13] 



The isobar degrees of freedom may of course also be considered explicitely in the solution of 
the nuclear many-body problem, i.e. one allows for many-body wave functions which contain 
isobar configurations. Such investigations have been made e.g. for the three-nucleon problem 
T9| and for nuclear matter [T3J. Presently, there is a renewed interest in the study of these 



isobar configurations in the nuclear wave function. One reason is the attempt to explore 
the isobar components in the nuclear wave function by means of different experiments. This 
includes the photoproduction of pions [[14] and the isobar current contribution to exclusive 
(e, e'NN) reactions. Another reason is the development of new models for the baryon-baryon 
interaction including isobar configurations ||l£)|-|T8|j . 

In this study we would like to investigate the predictions for the isobar components in the 
nuclear wave functions derived from various models for the baryon-baryon interaction. We 
are going to compare the isobar amplitudes calculated for the Argonne V28 potential (16 



with those derived from a recent update JT8| of Model I in table B.l of [2C|| , which we 
will denote as Bonn2ooo- Both of these models are based on the meson exchange picture 
of the baryon-baryon interaction. This means that e.g. the transition potential NN — > 
NA is dominated by the contribution from n exchange. While the Bonn 20 oo approach 
considers the complete relativistic structure of this term, the V28 model reduces this n 
exchange contribution to a local potential. This 7r exchange contribution is supplemented 
by the p exchange in the Bonn2ooo model. The V28 approach on the other hand adds a 
phenomenological contribution of short range. In [^T| it has been observed that the reduction 
of meson exchange interactions to the local approximation yields characteristic differences. 
Can such differences also be observed for the isobar contributions? 

A quite different approach has recently been developed by the Salamanca group [[17]] . They 
derive a NN interaction in the framework of the Chiral Quark Cluster (CQC) model. The 
problem of two interacting clusters (baryons) of quarks is solved by means of the resonating 
group method. The Pauli principle between the interacting quarks is an important source 
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for the short-range repulsion of the NN interaction f22| . At large distances the 7r exchange 
between the quarks in the two clusters evolves to the ir exchange between two baryons. At 
shorter distances, however, this non-local model for the baryon-baryon interaction might 
yield results that are quite different from a meson-exchange picture. This Salamanca poten- 
tial does not give such a perfect fit to the NN scattering phase shifts as the Bonn2ooo or the 
V28. For the 1 Sq and 3 Si — 3 D\ partial waves of the NN system, however, the agreement 
with the empirical data is rather good. 

We are going to compare correlated two-baryon wave functions in the deuteron, 16 O and 
infinite nuclear matter derived from these three models of the baryon-baryon interaction. 
While the deuteron wave functions are obtained from an exact diagonalization in momentum 
space, the wave functions for 16 and nuclear matter are calculated in the framework of the 
"exponential S" approach |2[], generalized to include isobar degrees of freedom and evaluated 
in momentum representation. 

After this introduction we discuss some differences between the transition potentials for V^28 
and Bonn 2 ooo m the next section. The wave functions of the deuteron derived from these 
three models are compared in section 3, while the correlation functions for 16 O and infinite 
nuclear matter are presented in section 4. A summary and conclusions are added in the final 
section. 



II. TRANSITION POTENTIALS IN DIFFERENT MODELS 

The key for an understanding of isobar components in the many-body wave function of 
nuclear systems is the transition potential, i.e. the matrix elements of the baryon-baryon 
interaction connecting two-nucleon states to A^A and AA configurations. A transition N <-> 
A implies the change of both spin and isospin from | to |. As isospin conservation has to be 
fulfilled at each vertex, such a transition requires the emission or absorption of an isovector 
meson, if one is considering a meson-exchange model for the baryon-baryon interaction. The 
most important isovector meson contributing to such transitions is of course the pion. The 
p(770) has a considerable effect, too, but mainly for large momentum transfers between the 
interacting baryons. 

In this section we would like to compare the one-pion-exchange (OPE) contribution as it 
is included in the Bonn 2 ooo interaction with the corresponding term in the Argonne V2& 
interaction. As an example we consider the transition potential AW «-> A A, in particular 
we will focus the attention to the central component as it shows up in 3 Si partial waves. 
This OPE transition amplitude has been formulated in a helicity representation in which 
the quantisation axis of the particle spins are the directions of the relative momenta q and 
q' for the initial and final state, respectively 



(q' AjVi Ajv 2 \V n \<\ Aa]Aa 2 ) = j^y 3 ^c(T) « - q,){q' v - q v ) i^q' - q) 

u(-q', AjvJ u"(-q, A A2 ) u(q', A. Nl ) u"(q, A Al ] 



q' - q) 2 + ml J (q' - q) 2 + m 2 + m A - m N 



c(T) is the isospin factor that remains from isospin matrix elements and corresponds to 
c(T) = — v2 in this channel. The energy denominator in this expression contains the 
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masses of the pion m^, the nucleoli and the tua- It deviates from a simple pion-exchange 
propagator, as it has been derived from time-dependent perturbation theory accounting for 
the mass difference between and m& |[24j| . The function F n represents the form factor 



for each ttNA vertex. The initial relative momentum q and the final relative momentum q' 
are chosen along the z axis, q = (0, 0, 0, q), and in the x-z-plane, q' = (0, q' sin 9, 0, q' cos 9), 
respectively, where 9 is the angle between q and q'. labels the helicity quantum numbers 
of the spinors. For the A isobars these spinors fulfil the Rarita-Schwinger equations and can 
be constructed explicitely by coupling a vector field (spin 1) and a Dirac field (spin |) 

-A 

Mq> A ) = £ c iM s e M> x ) s ), ( 2 ) 



A,s 



with 



t w i eA ' q a q e A-q , 

^ = ( — -> " e A , 1 , (3) 

'"A mAl v q + + m A) 



where the CAi are Clebsch-Gordan coefficients in the notation of |25f and w(p, s) is a Dirac 
spinor in spin state s. e+, eo and e_ are the circular polarisation vectors. 
From both the Dirac spinors for the nucleons and the Rarita-Schwinger spinors for the Deltas, 
the complex momentum structure is now removed by taking their value at q = q' = 0. It is 
straightforward to obtain the helicity Feynman graph structure (|l|) for all 2 x 2 x 4 x 4 = 64 
combinations of helicity projection numbers in this 'static' approximation. Most of them can 



be derived from symmetry relations |23|, and one only needs to evaluate 10 matrix elements 
explicitely. 

The helicity amplitudes are then projected onto states with definite total relative angular 
momentum J: 



(q 1 J A 'Ajv 1 Ajv 2 |V 7r |5 «/ A AaiAa 2 ) = 27ry^d(cos^)rf^ jA /(^)(q / Aiv 1 AAr 2 |y 7r |qAAiAA2)- 



(4) 



Here, d J KK ,(9) are the reduced rotation matrices with A = Aa x — Aa 2 and A' = Ajv x — Aat 2 . 
In the next step one performs the transformation into the basis of the partial wave states, 
in which orbital angular momentum L and spin S are coupled to J: 



{q Lj lVW Lj) - 2J + 1 ^a'CWj-a* 

AjVj A]v 2 Aa 1 Aa 2 



(g , J A 'A^A JVa |V'| ff J A AA 1 AA.>C7fi Ai| _ AAfl Ca A . (5) 



The transformation coefficients have been taken from [26]. 

In order to compare these matrix elements for the transition potential in the Bonn2ooo model 
with the corresponding one for the V28, we consider two approximations: 
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In evaluating the helicity amplitudes of eq.(|l]) the non-relativistic static limit is con- 
sidered, i.e. the Dirac spinors as well as the Rarita-Schwinger are considered for 
q = q' = 0. 

The energy denominator in eq. (0) is replaced by the usual pion propagator by ignoring 
the mass difference — m^. 



Using these approximations an analytical expression is obtained for 



'•VV.V|("|„:ic-i-i\ - ^ fNAir f 1 j(„^q\ ^_ m l 



( 2n) 9 mi J-i 



q' - q) 2 + ml 



^ 2 (q'-q) (6) 



that depends on the momentum transfer q' — q only. This means that the interaction 
is local. Transforming this expression for the central part of the transition potential into 
configuration space, we can identify a Yukawa term and a contact term, both multiplied by 
the form factor. 

This expression can be compared with the local transition potential as defined for the Ar- 
gonne V28 interaction, 

(r 3 S? N \V*\r 3 Sf A ) = ( f " A « m A e —L(\ - e - cr2 )fS'f JV |(S 1 S 2 )+(T 1 T 2 )t| 3 i S 1 AA ), (7) 

\ Viz J in~ r 

with Si and Ti denoting the transition operator for spin and isospin. Note that comparing 
this expression with the local approximation to the Bomi2ooo model in the contact term 
is removed and the form factor is replaced in V28 by a Gaussian cutoff to regularize the 
matrix elements for small interparticle distances. 

Matrix elements of the one-pion-exchange (OPE) contribution to the transition potential in 
the partial wave (q' 3 S^ N \V n \q 3 S AA ) for these different approaches are displayed in Fig. |l|. 
The relative momentum of the NN pair q' is fixed at 100 MeV and the matrix elements are 
presented as a function of the relative momentum for the AA state. 

The mass difference m& — in the pion propagator of ([!]) has a remarkable effect on the 
calculated transition potential. This can be seen in Fig. [I] from the comparison of the solid 
line (Bonn), which represents the complete OPE contribution in the Bonn 2 ooo potential, with 
the long dashed line (7r-range), which exhibits the results obtained after replacing the pion 
propagator of (fj) by the conventional n propagator. The inclusion of the mass difference 
ira — m^r leads to a quenching of the transition potential by about a factor two for small 
momenta, while the two curves approach each other at large values of ^aa- This means that 
the inclusion of the mass difference leads to a transition potential that is generally weaker 
and of shorter range, a feature which has already been observed e.g. in ||27|| . 
If one furthermore employs the non-relativistic limit, one obtains the local representation 
of the OPE of (|6]), which is given in Fig. [L] by the dashed line (local). The removal of 
the relativistic features yields a sizeable effect at larger momenta in particular. The dashed 
dotted line (contact rem.) results from the expression of (|GJ) if the contact term, the constant 
in the momentum representation, is ignored. The removal of this contact term has a very 
strong effect on the transition potential in this central partial wave with / = 0. Only after 
removing this contact term, we obtain an OPE component of the Bonn potential, which is 
essentially identical to the OPE contribution in the Argonne V28 model. 
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FIG. 1. Matrix elements for the one-pion-exchange contribution to the transition potential 
(q' 3 S | V 71- 1 g 3 S'^ A ) assuming various approximations as discussed in the text. Results are dis- 
played as a function of the momentum q in the AA state. The momentum of the NN state has 
been fixed to q = 100 MeV. All values have been multiplied with a common factor of 10 6 . 

The comparison of these various approximations in Fig. [I] makes it rather obvious that the 
different treatment of the OPE contribution to the NN — > AA transition potential leads to 
quite different matrix elements in the relativistic Bonn model, denned in momentum space, 
as compared to the local treatment in the Argonne V28. 

This OPE contribution is the most important ingredient to the transition potential in both 
models. The addition of the p exchange in the case of the Bonn 2 ooo model leads to minor 
although non-negligible modifications. This can be seen from the comparison of Fig. [1] with 
Fig. which shows matrix elements of the total transition potential. Besides the results for 
the Bonn2ooo and V28 models, this figure also shows the corresponding values for the quark 



model of the Salamanca group ||17|| . We see again the huge differences between the V28 and 
the Bonn 2 ooo model, which even lead to a difference in sign. From our discussion above we 
know that the main part of differences can be traced back to the different treatment of the 
OPE contribution. 



The Salamanca Chiral Quark Cluster (CQC) model, which we study as a third example fl7| 



cannot be considered as an entirely realistic potential because it does not fit the P-phase 
shifts for NN scattering with good accuracy; the phase shifts in the 1 <5'o and 3 Si — 3 Di 
channels are reproduced very well, however. In the CQC model, the long-range interaction 
is generated by OPE as well. In the core region, the potential is determined by gluon 
exchange and the Pauli principle. The latter requires total antisymmetry of the six-quark 
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FIG. 2. Matrix elements for the transition potential NN — > AA in the 3 Si central partial wave 
assuming the Bonn2oo(b the Argonne V28 and the Salamanca quark model approach. Further 
details see Fig. [p. 

wave function and this leads to a dependence of the baryon-baryon potential on the initial 
kinetic energy E in of two interacting nucleons. For the results displayed in Fig. ^| we assume 
a value E in = for this energy. 

The matrix elements for the transition potential of the Salamanca model exhibit values in 
between the two approaches discussed above. The shape is similar to the one of the Bonn 
potential, however, approaching the value of zero with increasing momenta much faster than 
the Bonn potential. This might be an indication that the Pauli effects in the quark model 
provide a much stronger cutoff at high momenta than the form factor used in the Bonn 
potential. As another example we compare results obtained for these three models also 
for the case (q' 3 Si N \V\q 7 D AA ) in Fig. [| The differences between Bonn 2 ooo an d V28 can 
again be traced back to the different treatment of the OPE. For this partial wave, however, 
one does not obtain any contact term contribution. This implies that the V28 is stronger 
at small momenta as it ignores the mass difference — in the pion propagator and 
weaker than Bomi2ooo at high momenta because of the non-localities included in the Bonn 
model. The matrix elements for the Salamanca model are similar to those of the Bonn 20 oo 
potential at small momenta but exhibit features of a stronger cut-off at high momenta. 
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FIG. 3. Matrix elements for the transition potential (q' 3 S^ N \V\q 7 D AA ) assuming the Bonn2ooo> 
the Argonne V28 and the Salamanca quark model approach. Further details see Fig. |l|. 

III. THE WAVE FUNCTIONS OF THE DEUTERON 



Solutions of the bound state two-particle problem have been obtained for the Argonne V28, 
the Bonn 2 ooo an d the Salamanca CQC potential. The nucleonic part of the deuteron wave 
function contains components 3 Si — 3 Di partial waves. If A degrees of freedom shall be 
taken into account explicitely, one has to extend the two coupled nucleonic channels by four 
AA partial waves, namely 3 S AA , 3 D AA , 7 D AA and 7 G AA . No A^A states can occur because 
the deuteron is an isospin T = state and the isospins of A" and A cannot couple to zero. 
The problem to determine the baryonic wave function in the 6 coupled channels has been 
solved in a spherical box with radius Rb ox = 20 fm. The spherical Bessel functions with the 
boundary condition 

ji(k u R box ) = 

can be used to construct a complete basis of orthonormalized states for the bound state wave 
functions within the spherical box [pB ]. The coupled channel Hamiltonian was calculated 
in this basis of momentum eigenstates of the box and diagonalised numerically. Up to 300 
discrete momenta were needed to get stable results for the binding energy. The results for 
energies and the wave functions obtained with the Argonne V28 and the Salamanca potential 



are in good agreement with the values given in |16j and ||17|| , respectively. 

All AW interactions, which we consider, are adjusted to fit the total energy of the deuteron, 

which can be written as a sum of the kinetic energy T tota i and potential energy V to tai with 



total 



NN 3 S 1 \T\NN 3 S 1 ) + (NN 3 D 1 \T\NN 3 D 



(AA|T|AA) 
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Arg. V28 


Bonn 20 oo 


Sal. CQC 


Arg. V14 


CD-Bonn 


Tfotal 
Vtotal 


23.75 
-25.97 


22.29 
-24.51 


17.82 
-20.04 


19.22 
-21.44 


15.48 
-17.70 


rpN 
1 


10.32 
8.95 

A AO 

4.48 


10.25 
5.35 
6.69 


10.84 
5.13 
1.85 


10.54 
8.68 


9.79 
5.69 


V SS 

V DD 
i/NN 
V SD 
yNA 

yAA 


7.08 
5.86 
-29.22 
-10.40 
0.71 


1.52 
1.80 
-14.38 
-13.44 


-5.67 
0.70 
-11.64 
-3.16 
-0.28 


-1.83 
1.99 
-21.60 


-4.77 
1.34 
-14.27 



TABLE I. Contributions to the kinetic and potential energy of the deuteron originating from 
the different parts of the wave function. The kinetic energy Ttotai is the sum of the kinetic energies 
originating from the iViV wave function in 3 Si (Tg) and 3 D\ (Tp) partial waves plus the kinetic 
energy (T A ) from the AA components. The term T A also accounts for the NA mass difference. 
The potential energy contains contributions from the various parts of the NN —* NN potential 
(V NN ), the NN -> AA terms (V NA ) and the AA -> AA terms (V AA ). Results are given for 
the 3 models with isobar configurations. For comparison we also show results from pure nucleonic 
interaction models V14 and CDBonn. All entries are given in MeV. 



= Tg + T% + T A . (8) 

Note that the term T A sums up the kinetic energy from all partial wave components in 
the AA wave function and includes the contribution from the N — A mass difference. The 
potential energy can be split into 



V to tai = (NN 3 S 1 \V\NN 3 S 1 ) + (NN 3 D 1 \V\NN 3 D 1 ) + 

2 (NN 3 Si\V\NN 3 D 1 ) + 2 (NN\V\AA) + (AA|V|AA) 
= V S N S N + vSS + V S N D N + V NA + V AA . (9) 

The contributions of these various terms to the energy of the deuteron are listed in table |, 
while the norm of the various partial wave contributions to the wave function are presented 
in table |I| for the three interaction models with inclusion of isobars. For a comparison 
table | also shows the result for two NN interaction models without the explicit treatment 



of isobar configurations: the Argonne potential model V14 and the charge-dependent 
Bonn potential 0. 

Although the sum of these various energy contributions yields the same energy of the 
deuteron for all the interaction models, there are remarkable differences in the individ- 
ual terms. The inclusion of isobar components in the wave function yields a contribution 
to the binding energy of -5.21 MeV, -6.75 MeV and -1.59 MeV for the V28, the Bonn 20 oo 
and the Salamanca interaction, respectively. This means that V28 and Bonn 2 ooo predict an 



9 







V 28 



Bonn 20 oo 



CQC 



3 gNN 
3 D NN 

3 S AA 

3 D AA 
7 D AA 
7 qAA 

Total AA 



93.341 
6.133 
0.043 
0.020 
0.417 
0.045 
0.524 



94.685 

4.705 

0.223 

0.022 

0.361 

0.005 

0.611 



95.199 
4.567 
0.107 
0.004 
0.124 
0.006 
0.241 



TABLE II. Probability for D-state and AA-components in the deuteron for different models of 
the baryon - baryon interaction. All entries in percent. 



unbound deuteron, if the AA components in the wave function would be suppressed, while 
the Salamanca interaction leads to a much weaker contribution of the isobar configurations. 
These differences also show up in the norm of the various partial wave components in the 
deuteron wave function listed in table ||. Comparing these occupation probabilities one finds 
that all interactions predict the largest contribution to the AA probability in the 7 Di partial 
wave. The largest probability for this partial wave is obtained for the V28 potential, while 
Bonn2ooo and the Salamanca interaction in particular yield much smaller probabilities in this 
channel. The origin of these differences can be deduced from Fig. |], which exhibits the AA 
wave functions. The Argonne potential leads to a wave function in the 7 Di, A A channel, 
which is of longer range than those derived from the other two interactions. The main 
reason for this feature is the fact that the OPE contribution to the NN — > AA transition 
potential neglects the N — A mass difference in the pion propagator (see ([I])), which leads 
to a transition potential of longer range. Comparing the wave functions in this channel 
derived from the Bomi2ooo and the Salamanca interaction, one observes that the latter is 
strongly suppressed at short relative distances as compared to the former, which supports 
the argument presented already in the comparison of the matrix elements of the interaction, 
that the quark model leads to a much stronger effective cut-off at short distances than the 
cut-offs used in the meson exchange model. 

The most significant difference in the predictions for the isobar components of the wave 
function can be seen in the 3 Si, A A channel. While the Bonn and Salamanca models yield 
occupation probabilities, which are of similar size as those for the 7 -D 1; A A channel, the 
probability for this channel derived from the V28 potential is much weaker. A qualitative 
difference can also be observed from the inspection of the wave function in this channel 
(see Fig. [|). The wave function obtained for the Bonn and Salamanca potential show a 
maximal amplitude for r — > 0, while the corresponding wave function deduced from the V28 
potential is suppressed at small relative distances. These differences in the wave function 
reflect the fact that the contact term in the OPE contribution to the transition potential 
(see (^|) and discussion there) is removed in the local representation of the OPE in V28, 
while the non-local representation of the OPE in the other models keeps a strong short 
range component. 
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FIG. 4. Isobar components in the deuteron wave function using the three different models for 
the baryon-baryon interaction. 

The comparison of the AA components in the deuteron wave function therefore reflects the 
main differences between the local approximation of the V28 and the non-local approaches 
of the Bonn and Salamanca model: The V28 contains a OPE component of longer range 
and suppresses the short range components much stronger than the other two models. Fur- 
thermore, the quark model approach leads to a stronger reduction than the typical cut-off 
that is used in the meson exchange models of the Bonn group. This leads to much weaker 
AA components in the deuteron wave function for the Salamanca model. 
It is worth noting that the contributions of Vg^ to the potential energy of the deuteron 
(see table |) are significantly larger for the local Argonne potentials V14 and V28, than for 
the various versions of the Bonn potential and the Salamanca potential. This observation 
suggests that the local representation of the tensor components of the NN interaction, in 
particular the contribution originating from 7r-exchange yields matrix elements which are 
significantly larger than those evaluated within the non-local approaches (see also 




IV. CORRELATIONS IN NUCLEAR MATTER AND ie O 

The wave functions of many-body systems are studied in the framework of the coupled- 
cluster theory |J. In the coupled-cluster approach one starts assuming an appropriate 
Slater determinant $ and writes the exact eigenstate $ for the v4-particle system as 

* = e 5 $, (10) 

with S an operator of the form 

s = £s n , (li) 

n=l 
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where S n is an n-particle operator which can be written for the case of n = 2 

S i = \ (PiP2\S 2 \viv 2 ) a ] pi al 2 a V2 a Ul . (12) 

In this equation aj,. stand for fermion creation operators in states which are unoccupied in 
<3>, while a Vi represent annihilation operators for the nucleon single-particle states which are 
occupied in the Slater determinant $. Note that the aj,. may also represent the creation 
of A isobar states. Therefore the Si amplitudes describe two-particle two-hole excitations 
relative to $ but also A A and AA excitations. 

For our present investigation of nuclear matter we will assume <3> to be the Slater determi- 
nant defined in terms of plane waves, occupying all states with momenta up to the Fermi 
momentum kp = 1.36 fm^ 1 . As an example for a finite nucleus, we will also consider 16 0. In 
this case we will assume $ to be the Slater determinant, defined in terms of harmonic oscil- 
lator states (hu = 14 MeV) with nucleons occupying the states of Os and the Op shell. If one 
assumes that these single-particle states represent an optimal single-particle basis, i.e. the 
amplitudes Si in (|ll]) vanish, and ignores the contributions of linked n-particle correlations 
with n > 3 (S n = for n > 3), one obtains integral equations for amplitudes 

(6162 [k(lS)j] KLJtISzKv^Jt) , (13) 



which are solved in momentum space |3(J. In this representation &1&2 stand for AA, A A 



and A A states, k(lS)j denote the momentum, spin and orbital quantum numbers for the 
partial wave basis of the relative motion of the two baryons. K and L represent the center 
of mass state and J and r refers to the total angular momentum and isospin of the pair of 
baryons. 

The 5*2 amplitudes can be considered as correlation functions describing the difference be- 
tween the correlated and uncorrelated wave function of two particles in the nuclear medium. 
As the uncorrelated Slater determinant $ in (|i~0| ) does not contain any isobar components, 
the £2 amplitudes can be interpreted directly as the A A and AA component of the two- 
particle wave function, if 6162 in fljjD refer to A A and A A states. 

As typical examples for these correlation functions in nuclear matter, the AA components 
for the relative wave function of two baryons in a state with isospin r = and angular 
momentum j — 1 are given in Fig. |5|. These components correspond to the components 
of the deuteron wave function in the same partial waves displayed in Fig. In fact, these 
isobar components of the wave function for a pair of baryons in nuclear matter is very 
similar to the deuteron wave function if one compares the partial wave with I = 2. The AA 
wave functions exhibit a tail of longer range in the case of the Argonne V28 interaction, 
which reflects the longer range of the AA — > AA transition potential for this interaction as 
compared to the other two approaches. The Salamanca CQC approach yields amplitudes 
for these D-waves which are considerably smaller. All these amplitudes are enhanced as 
compared to the deuteron wave function, which reflects the larger density of the nuclear 
matter system. 

The situation is a little bit different if one compares the 3 Si, A A component of the nuclear 
matter wave function with the corresponding part of the deuteron wave function. The results 
obtained for the V28 and the Bonn potential yield wave functions of rather similar shape 



12 



Nuclear Matter, x=0, AA 




FIG. 5. Isobar components in the relative wave functions of two particles in nuclear matter with 
Isospin t = and J = 1. 



in the deuteron and in nuclear matter also for this channel. The wave functions derived 
from the Salamanca CQC interaction are rather different in nuclear matter as compared 
to the deuteron. At short distances they even have a different sign. This reflects a large 
non-locality or momentum dependence of the short range component for the NN — > AA 
transition potential derived in the CQC model. 

In the nuclear many body systems one also observes bound states of two baryons in states 
with isospin r = 1. As an example of such configurations we discuss NA and AA compo- 
nents of two-particle wave functions in 16 as displayed in Fig. [| The comparison of the 
results obtained for the different interaction models leads to observations which are rather 
similar to those discussed for the states with r = above. The difference in the range 
of the transition potentials iViV — > NA for V28 and the other two interactions is not as 
significant as in the case of the NN — > AA (compare the discussion of the propagator in 
(H)), which leads to smaller differences in the tail of the 5 Do NA wave function than in the 
corresponding AA state. A very significant model dependence is obtained in the 1 S'o, AA 
wave function. 

The isobar-isobar relative wave functions derived for nuclear matter and finite nuclei are 
rather similar. Therefore we display only one of these examples for each channel. This 
demonstrates that isobar admixtures correspond to correlations in the many-body wave 
function which are of short range. Therefore they are not very sensitive to surface effects 
in finite nuclei and a local density approximation should be appropriate to consider isobar 
effects in finite nuclei. 
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O, Isobar components, x=l 
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FIG. 6. Isobar components in the relative wave functions of two particles in 16 with Isospin 
r = 1 and J = 0. The uncorrected state corresponds to two nucleons in the Osi/2 state. 



V. CONCLUSIONS 

Three different models for the baryon-baryon interaction, which fit NN scattering data and 
explicitely account for isobar degrees of freedom have been considered. These are the local 
Argonne V28 potential [16], the Bonn2ooo [18 interaction model based on the relativistic 



meson exchange model and an interaction based on the Chiral Quark Cluster (CQC) model 
which has recently been developed by the Salamanca group [17]]. The isobar components 
in the wave function of the deuteron and nuclear many-body systems including nuclear 
matter and 16 O as an example for a finite nucleus are evaluated. Significant differences 
are observed in the predictions from these models. These differences can be related to the 
assumptions made in determining the transition potential between iViV, NA and AA states. 
The V28 interaction model yields isobar wave functions of longer range than the other two. 
This can be traced back to the approximations which are made in reducing the 7r exchange 
contribution of the transition potentials to a local form. The most significant differences are 
observed in partial waves with I = at small distances. A main source of the discrepancies 
in these partial waves is related to a removal of a contact term in the local interaction 
model. The Bonn 2 ooo an d Salamanca CQC model in particular exhibit important non-local 
features for the short range part of the interaction model. The Salamanca model predicts 
rather small isobar components in the nuclear many-body wave function. These different 
predictions for the isobar components in the many-body wave function could be very useful 
in distinguishing between different interaction models. 
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